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<n : 

O ' Abstract 

<d : 

Q . We construct universal parton evolution equation that produces space- and time-like anoma- 

I lous dimensions for the maximally super-symmetric M = 4 Yang-Mills field theory model, and 

find that its kernel satisfies the Gribov-Lipatov reciprocity relation in three loops. Given a sim- 
ple structure of the evolution kernel, this should help to generate the major part of multi-loop 
contributions to QCD anomalous dimensions, due to classical soft gluon radiation effects. 

00 1 

1 Introduction 

\Q \ Super-symmetric Yang-Mills theories, which neighbour QCD in their physical content but are 

easier to handle, may play an important role in shedding light upon multi-loop QCD results. 

4h ■ In this work we employ the J\f = 4 super-symmetric Yang-Mills (M = 4 SYM) field theory in 

order to better understand some important features of the QCD parton evolution. 

<D ! 

Soft gluons and the "inheritance" idea. Permanent fight for increasing precision of 
QCD predictions necessitates higher order perturbative calculations. Such calculations, and 
markedly the recently completed three loop analysis of DIS anomalous dimensions [1], yield 
lengthy expressions reach in hidden physical content that one has to grasp. A higher loop 
expression necessarily contains certain structures inherited from the preceding orders of the 
PT expansion. This is the case, for instance, for the renormalization effects that build up the 
physical coupling a p h determining in all orders the intensity of soft gluon emission given by 
the Born level expression [2,3]. 

We expect this "inheritance" feature to go deeper and to include, in particular, real and 
virtual radiative effects due to multiple soft gluons. The basis for such expectation is the 
celebrated, virtually unknown, Low-Burnett-Kroll theorem [4] which fully applies to QCD as 
well as to SYM quantum field theories. It tell us that soft gluon radiation has in fact classical 
nature. It is independent of the quantum state of the radiating parton system, and does not 
change it. Therefore it should be possible to find a procedure for iteratively "dressing" an 
underlying lower order hard parton splitting process by soft gluons. The realisation of this 
programme would allow one to generate, by simple means, the corrections that the anomalous 
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dimension receives from soft gluons in all orders of the perturbative expansion. There is a 
good reason to believe that soft gluon radiation effects are responsible for the major part of 
high order contributions to the QCD anomalous dimension. At the two loop level this has been 
explicitly demonstrated in a particular example of the time-like heavy quark fragmentation 
function in [2]. 

To verify this expectation at the three loop level we use as a testing ground the M = 4 
SYM model. Since soft radiation is universal, the results of the present study are directly 
applicable to QCD (modulo adjustment of colour factors). Indeed, the diagonal one loop 
anomalous dimensions (q — > q, g — > g) contain the universal classical part of the gluon emission 
probability, 



— - C ^ 

dx P 7T 



^— - + (1 - x)g p (x) 



(1.1) 



Being proportional to the "colour charge squared" C p , it is otherwise insensitive to the nature 
of the radiating particle p. Quantum effects encoded in g p obviously depend on p (g±/2 = \ 
for a fermion, g\ = x + x _1 for a vector emitter, etc.). Importantly, as a consequence of the 
LBK theorem, the genuine quantum contribution to (jl.ip is down, relative to the universal 
classical piece, by two powers in the (1 — x) counting. In the standard approach, in higher 
orders classical and quantum effects get mixed by the Q 2 evolution, resulting in cumbersome 
expressions for anomalous dimensions. 

In the Af = 4 SYM theory, g p = in (|l.ip at one loop: quantum effects due to boson and 
fermion fields compensate each other, leaving us with a pure "classical physics" at the Born 
level, and thus provide a perfect ground for testing the "inheritance" idea. The very idea of 
extraction of the third loop contribution to the Jsf = 4 SYM anomalous dimension [5] from 
the corresponding QCD result [1] was based on the observation that the two-loop expression 
for the anomalous dimension 7sym is built of the functions characterised by the same level of 
transcedentality [6]. In fact, such functions (with transcedentality r = 2k — 1 in the k th order 
of perturbation theory) are nothing but finite "leftovers" of the real-virtual cancellation of 
infrared gluon divergences. Bearing this in mind, the study of higher orders in the "universal 
anomalous dimension" of the M = 4 super-multiplet of twist-two operators [5] should help us 
to understand the structure of the iteration of pure classical gluon emissions in QCD. 



Analytic link between parton distribution and fragmentation functions. It 

is important to stress that apart from being large, soft gluon effects are also responsible for 
generating enhanced ln n (l — x) corrections in the quasi-elastic kinematics, l-i < 1, both for 
parton distribution functions in space-like process, xb — > 1, and fragmentation functions in 
time-like parton cascades, xp —* 1- Such corrections translate into IniV terms in the asymp- 
totics of anomalous dimensions at large N. In one loop [7-9], where 7(JV) — C p ^-ip(N) oc In TV, 
the task of analytic continuation was solved in [10] where it was proposed to treat certain log- 
arithms in arithmetic sense, ln(l — x) => In 1 1 — x\; in the QCD context this idea had been 
implemented in [9]. However, in higher orders the power n of the (lniV) n enhancement in 
the anomalous dimension increases, thus impairing the analytic link between DIS and e + e~ 
annihilation channels. For a thorough discussion of the problem at the two-loop level see [11] 
and references therein. 

The Af=4 SYM framework is perfectly suited for shedding light on the problem of finding 
a way to analytically relate these cross-channel processes. For the benefit of the reader, let us 



2 



recall this problem in more detail. 

In the DIS case the large virtual momentum q transferred from the incident lepton to 
the target nucleon with momentum P is space-like, q 2 < 0. Inelasticity of the process is 
conveniently characterised by the Bjorken variable xb = —Q 2 /2(Pq)- Inclusive fragmentation 
of an e + e~ pair with total momentum q (large positive invariant mass squared q 2 ) into a final 
state hadron with momentum P is characterised by the Feynman variable xf = 2{Pq)/q 2 
(hadron energy fraction in the e + e~ cms.). The fact that Bjorken and Feynman variables are 
indicated by the same letter is certainly not accidental. In both channels < x < 1 though 
these variables are actually mutually reciprocal, xf 1/xb, rather than identical. One x 
becomes the inverse of the other after the crossing operation P^ — > — P^: 

~q 2 2(Pg) 
x B = x F = —. (1.2) 

Apart from the difference in the hadron momentum P belonging to the initial state in the DIS 
and final state in the e + e~ case, the Feynman diagrams for the two processes are basically the 
same. "Mass singularities" that emerge when additional momenta in the Feynman diagrams 
become collinear to P are therefore also the same. That is why in the two processes a similar 
parton interpretation emerges in terms of QCD evolution equations. 

The common structure of Feynman diagrams and the reciprocity (11.20 are the origin of the 
close relation between deep inelastic structure functions and e + e~ annihilation inclusive frag- 
mentation functions. In particular, the space- and time-like evolution anomalous dimensions 
turn out to be related. 



Drell-Levy-Yan relation. Back in 1969, Drell, Levy and Yan have proposed that 
the two channels can be linked by analytical continuation [12]. If one takes the space-like 
splitting function (x) — the Mellin transform of the anomalous dimension 7^) (N) — and 
analytically continues it from the physical region x = xb < 1 to x > 1, such a procedure would 
yield the time-like splitting function with xf = x < 1: 

7 (T)( S ) = -x" 1 7 {S) (x- 1 ) . (1.3) 

Gribov-Lipatov reciprocity. Gribov and Lipatov in their first systematic study of DIS 
and e + e~ annihilation processes in the QFT framework [7], suggested that the common dy- 
namics in the two channels results simply in the identity 

^)(x) = 7( s )(x). (1.4a) 

Taken together with the DLY relation fjl .31) this leads to the relation known as the GL reci- 
procity: 

7(x) = -17 (x~ l ) . (l-4b) 

In the Mellin moment space, 

7 (iV) = M[7(x)]= / -x N j(x), (1.5) 
Jo x 

(jl.4p translates into the symmetry iV — > —(N + 1); in other words, in order to satisfy the GL 
reciprocity, the anomalous dimension "f(N) has to be a function of a single variable N(N+1) — 
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the Casimir operator of the collinear 50(1, 2) subgroup of the 50(2,4) conformal group [13]. 
A group symmetry nature of the GL reciprocity was advocated already in the early 70's when 
the issue was first raised [14]. Indeed, since x is a light-cone momentum fraction, y = Inx 
represents parton (pseudo)rapidity: ko = m±coshy, k z = m^sinhy. Therefore in the Mellin 
transformation, 



J ^L X N = I dye yN ^ j f 



an analogy between the Lorenz boost parameter y and an imaginary rotation angle allows one 
to look upon the conjugate variable N as "angular momentum", with J 2 = N(N + 1) the 
corresponding Casimir operator. 

True in one loop (leading logarithmic approximation), see e.g. (jl.ip . the relations (jl.4p are 
known to break already in the second loop [15, 16]. 

Reciprocity respecting evolution equation. This letter presents a development of 
the project that aims to reorganise the QFT parton evolution picture [17] in such a way as to 
preserve a close relation between space- and time-like parton dynamics. In [2, 18] a "reciprocity 
respecting" equation (RRE) for QCD anomalous dimensions has been suggested, based on a 
notion of the universal "evolution kernel" P(x,a), identical for two channels, that satisfies the 
Gribov-Lipatov relation in all orders: 

d T DW(N,Q 2 )= 7a (N,a)D( T / s \N,Q 2 ) = C - z N P(z,a) (N, z°Q 2 ) . (1.6) 

Jo z 

Here r = InQ 2 , while a = +1 for e + e~ annihilation parton fragmentation functions Z?( T ) 
(time-like evolution), and a = — 1 for DIS parton distributions (space-like). Using the 
Taylor expansion trick we obtain the formal solution of (|1 .6|) : 



d T D^ s \N,Q 2 ) = f 1 -z N P(z,a)z^D^/ s \N,Q 2 ). (1.7) 
Jo z 

At the three-loop level, this equation has been used to analyse the x — > 1 behaviour of the 
splitting functions j( T / s \x) in [18]. There the all-order relation C = —a A 2 between the 
magnitudes of the A(l — x)~ l and Cln(l — x) singularities was established, observed by Moch, 
Vermaseren and Vogt for quark and gluon space-like anomalous dimensions [1] . The prediction 
of (|1.7p for the time-like case was verified at the third loop by explicit calculation of the non- 
singlet anomalous dimension in [19]. 

The analysis of next subleading contributions in the x — > 1 limit attempted in [18] was 
complicated by the presence of scheme dependent contributions proportional to the (5 function^, 
starting from the constant term j(x) = 0(1), which contribution also had to be under full 
"inheritance" control, according to the LBK wisdom. It is indeed in the N = 4 SYM model, 
having (3(a) = 0. 

Here the RRE (jl.7p readily reduces to a pair of complementary relations which hold in all 
orders and for arbitrary iV: 

7 (T/S) (iV ) = p(AT + ( 7-7( T / s )(iV)), (1.8a) 
V(N) = 7 (T/s) (iV -a-V(N)). (1.8b) 



: We thank G. Korchemsky for communicating us an elegant solution of the problem encountered in [18]. 
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Expanding (jl.8ap results in 



7 = 7 ( T / S )[a] = V + a-fV+^ 2 V + a^'v + 0(a 5 ), (1.9) 

where dots indicate derivatives with respect to N and V = V(N). Since 7 = 0(a), this equation 
gives iteratively a tower of inherited higher order contributions to the anomalous dimensions 
in both channels. Solving (jl.9p perturbatively gives 



i( 7 (T) +7 (S)) _ V + R, R= \vV 2 + \V 2 f\ +0(a 5 ), (1.10a) 
i( 7 (T)_ 7 (S)) _ v = VV+[^ + lV 2 VV + \V z v]+0{a & ) . (1.10b) 



The difference of space- and time-like anomalous dimensions V in (|1.10b|) does not have a 
linear term in the evolution kernel V; therefore, knowing V at order k, one predicts V at 
the order k + 1. Moreover, since the induced terms in (jl.lOp contain powers of V, the lowest 
order evolution kernel V\ oc lnJV generates a tree of logarithmically enhanced terms, up to 
a k ln k N/N k+1 in order k of the perturbative expansion. 

If the large-iV singularity of the evolution kernel V turns out to be stable under pertur- 
bation, given V at order k, (|1.9j) would predict all logarithmically enhanced contributions 
f m (N) ■ ln m N, 2 < m < k + 1 at the next order k + 1, where f m (N) = 0(N- k - 1 ) are known 
functions that are regular at N = 00. 

We translate the space-like N = 4 SYM anomalous dimension recently calculated by 
Kotikov, Lipatov, Onishchenko and Velizhanin in three loops [5] into the corresponding evolu- 
tion kernel V of (II. 8D . We find that Vz is given by a rather transparent expression, much more 
compact than 73 , exhibits only a single power of the In iV enhancement at large N and satisfies 
the GL reciprocity (|1.4p . Grasping internal logic of the second and third order contributions 
to V should help to gain better understanding of the structure of multi-loop QCD anomalous 
dimensions. 

The paper is organised as follows. In section [2] we present the evolution kernel V in three 
loops. In section [3] we expose in detail essential properties of the answer, including important 
kinematical limits. In section U] we present our conclusions and discuss perspectives for future 
studies. 



2 Answer 

We construct the RR evolution kernel, 



C A a 



k=l 



k=l 



7T 



(2.11) 



employing the space-like anomalous dimension, 7 = 7 ( - s ^, for the A/" =4 SYM model. Prom the 
perturbative expansion of (jl.lOp . 



7i = 'Pi, 72 = P2 ~ V1P1 



73 = Pi 



V1V2 + V1V2 



+ 



V{?\ + \V\V\ 



(2.12) 
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Within our normalisation convention (|2.1ip . one has 7^ 



7, 



/4 fc , with 7 un i the anomalous 



dimension in the original notation of [5] . The latter is given in terms of harmonic sum functions 
S± a (N) and S_ P)0) ...(iV), see [20,21] and AppendixfAj It reads 



71 = --Si; 

72 = 5S3 + S1S2 + ( 2 5-3 + S1S-2 — S-2,1) ; 

73 = — — [S1S3 + 5 ^2 53 + 5l5 2 + fiSiSi] 

-5 : [45_ 4 + ±S* 2 + 25 2 5_ 2 - 65_ 3 ,i - 55_ 2 , 2 + 85_ 2 ,i,i] 
— ( 2 5 2 + 35^ )5_3 — S3S-2 + (5 2 + 2 < S' 1 )5_ 2i i + 125_2,i,i,i 

—6(5-3,1,1 + 5-2,1,2 + 5_2,2,l) + 3(5_4,i + 5-3,2 + 5-2,3) _ § 



(2.13a) 
(2.13b) 

(2.13c) 



In order to arrive at a compact expression for the third loop evolution kernel V3 we introduce 
the characteristic function <p, 



dip 



dz z — 1 



<p(z) 



z z + r 

which is invariant under the transformation z 



In 



In cosh 2 (i In z 



4z 1 2 

z , construct integrals 



(2.14a) 



$t(x) 



1 



r(r) 



dz 



(<p(x) - <p(x)) 



T-l 



(2.14b) 



(where subscript r marks transcedentality of the function), and define Mellin transforms 



Y. m (N) = (-1) N M 



l + x 



m— 1 



(*) 



(2.14c) 



Functions Y- m {N) are given by linear combinations 

Y-3 = S-3 - 2 5i,_ 2 ; y_ 4 = 5_4 - 2(V3 + 5 2 ,_ 2 ) + 4 Si,i,_ 2 ; 

5-5 — 2(53,_ 2 + 5 2 ,-3 + 5i,-4) +4 (52,1,-2 + 5i,2,-2 + 5l,l,-3) 



85i 1 



1,-2) 



where /mi crowns a sum with subtracted N = 00 value: 

5 3 (AT) = S s (iV) - S s (oo). 



Here 5^,_ p (A^) are complementary harmonic sums, 



5 ai , g2 „„, aw (iV) = (-1)" £ 



(sgn ai 



=AT+1 



A; 



lail 



E 

fc 2 =fci+i 



(sgn a 2 ) fcl 



k 



12! 



• E 



(sgn a n )* 



(2.14d) 



(2.15) 



(2.16) 



whose relation to standard multi-index harmonic sums is displayed in Appendix [A] The 
functions (I2.14dj) are built of harmonic sums with one, and the rightmost, negative index 
(therefore the notation K_ m ). Internal logic of the constructs (|2.14dj) is transparent and makes 
generalisation straightforward. 

In these terms, the J\f = 4 SYM evolution kernel, in first three orders of the perturbative 
expansion in the physical coupling, 



«ph 



(1 - \Q 2 a + %CW + ...), 



(2.17) 
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takes a compact form 



v 2 



- Si; 

2 S3 



K-3 + B 2 



2^-3-1- -D2, 

+ |y_ 5 + B 3 + c 2 • ^3 



2 

+ 5i 



(2.18a) 
(2.18b) 

(2.18c) 



Here B 2 = |C 3 , 5 3 



"5C2C3 



|C5- 



3 Discussion 



Let us list some characteristic features of the RR evolution kernel (|2.18j) . 

The evolution kernel (|2.18p (as well as anomalous dimensions) contains harmonic sums 
with positive and negative indices. These contributions have different origin and play markedly 
different roles in extreme kinematical limits. Harmonic sums bearing a negative index (indices) 
oscillate with N, see dAlc]) and (|2.14c|) : 

V(N) = V pos (N) + (-1) N -V ncg (N). (3.19) 

The presence of the oscillating factor (—1)^ calls for introducing two separate functions, ana- 
lytically continued from even and odd moments N. Indeed, in QCD one has two independent 
non-singlet anomalous dimensions, replacing (|3.19p by 

7 £° = P qg ± vm- ( 3 - 2 °) 

The two anomalous dimensions have different signature with respect to s <-> u (x —> —x) 
crossing. In J\f = 4 SYM "quarks" are Majorana fermions, making q and q physically indis- 
tinguishable. In the light-cone formulation of AT = 4 SYM [22,23], in the study of anomalous 
dimensions [24] the negative signature does not appear [25], since operators built of the unique 
superfield of the theory have vanishing odd moments. Therefore, in the M = 4 SYM context 
the only possible continuation to non-integer moments consists of dropping the factor ( — 1)^, 
as it has been done in [5], resulting in 7 un i = 7 < - + - ) . We should keep in mind, however, that 
both signatures are present in the QCD case, where the maximal transcedentality structures of 
(|2.18p are also present and are responsible for a significant part of higher order contributions. 



3.1 Inheritance 

Compared to the standard anomalous dimension (|2.13p . the presence of harmonic sums with 
positive indices in the RR evolution kernel (|2.18p is minimal: 

72 = 5S3 + SiS 2 + • • • , V 2 = \Sz H ; 

73 = -\S b - [SfS 3 + ±S 2 S 3 + S1S2 2 + §Si S 4 ] +■■■ , V 3 = -±S 5 + • • • , 

where ■ ■ ■ mark contributions of negative index sums. Thus, all non-linear combinations of 
Si in the anomalous dimension has been generated by the RRE and are therefore physically 
"inherited" from the first loop. 
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3.2 Negative index sums: relation to Feynman graph disconti- 
nuities 



Sums with one negative index contribute to the anomalous dimension starting from the second 
loop. Their appearance can be traced to non-planar diagrams that can be cut both in s and 
u, and thus have non-zero Mandelstam double spectral function, p su ^ 0. 

The diagrams for non-singlet QCD quark evolution shown in Fig. [T] can be related by 
the crossing transformation s <-» u ~ — s which, in the language of the Bjorken x variable, 
formally corresponds to the reflection x — » —x. It is therefore only natural that the second 




(a) (b) 



Figure 1: On the origin of negative index harmonic sums 

loop non-singlet — > q~i anomalous dimension described by the graph (b), is proportional, in 
x representation, to the first loop % — > qi splitting function evaluated at the reflected point, 



Pqq{x) = Op h (\C A - C F ) Pqq(-X) ■ $ 2 (», Pqq(x) 



1 + x 2 
2(1 -x) 



(3.22) 



The proportionality factor <&2{x) is a functiorH of transcedentality 2, defined in (|2.14bj) . It is 
important to stress that this factor appears, invariably, in second loop anomalous dimensions, 
multiplying the full first loop splitting functions pba(-x) in all parton transitions, B — > A, 
including non-diagonal ones, see [16]. Using three-loop results of [1], it is straightforward to 
verify that in the non-singlet QCD anomalous dimensions 7^ every harmonic sum (max- 
imal transcedentality or not) bearing a negative index is proportional to 2Cf — Ca (colour 
suppressed), that is has non-planar origin. 

Let us compare contributions of positive and negative index sums. For the second loop 
evolution kernel we have 



M 



2 S3 



M 



^3-(-l) 



N 



\\n 2 x)=Vl°\x), (3.23a) 
•* 2 (a0) = V^(x), (3.23b) 



with Vi(x) = x/(l — x) the one loop splitting function, the inverse Mellin image of V\{N) = 
—S\{N). The fact that these two terms originate from different cuts of the same diagram can 



2 called 52(0;) in the original papers [15], [16] 



s 



be seen from internal structure of corresponding contributions: 



[i 



[l+x] 



1 dz [1 + z] 

z 1 + z 

1 dz [1 - z] 

T i + z 



lnz, 
lnz. 



(3.24a) 
(3.24b) 



where we have used (I2.14ap to transform the integral 



dz ( (1 + a;) 2 z 



= / — In 



(1 +z) 2 x 



1 dz z - 1 , 
inz. 

x 2z z + 1 



Higher order positive index contributions, £>2k-i in I12.18j) a.b. with k the order of the pertur- 
bative expansion, have a simple origin and a simple structure: 



S 5 (N) 



1 d 2 



2! diV 2 
1 d 4 



4! dN A 



Si(N), 
Si(N), 



M 



M 



1 x In 2 x 
2! x - 1 ' 
1 x In 4 x 



4! x - 1 



(3.25a) 
(3.25b) 



Given that the negative and positive index contributions are related by the s <-> u crossing, 
and guided by the second loop example, a possibility to generate the negative index piece from 
its positive counterpart does not look implausible. 



3.3 Gribov— Lipatov reciprocity 

The Gribov-Lipatov reciprocity is respected by each piece entering the evolution kernel. In- 
deed, it is easy to see from (|A.1|) that the inverse Mellin images of harmonic sums S a with odd 
a, as well as negative index sums S-b with b even, satisfy the relation 

F{x) = -xF{x- 1 ) (3.26) 

(For the distribution M -1 [Si] this relation holds for x ^ 1.) The same is true for functions 
f_ m , see (EUD, 

•*m-l(aO, (3.27a) 



ivr 1 



-1) N Y_ 



i 



which transform according to (|3.26[) . due to 

Qrix- 1 ) = -$ r (x). (3.27b) 

Finally, an observation that (|3.26p is stable under convolution of two RR functions, concludes 
the proof that every item in (|2.18p satisfies the GL reciprocity relation (|1.4|) . 

In the moment space, as we have discussed above, one expects symmetry with respect to 
the substitution 

V{N) = V(-N - 1), (3.28) 

and thus the dependence on the combination iV(iV + 1). Indeed, in QCD, for example, the 
second loop non-singlet positive signature evolution kernel reads [2] 



Vt' + \x) = {\C A - C F )[p qq (x) ■ \\n 2 x + p qq (-x) ■ $ 2 | 



x 



w . , ( 3 - 29 ) 



l 



0(i(l+s)lnx + (l-s)) +\C A {1- 
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It contains "algebraic" pieces whose dependence on the "Casimir" is transparent, 

M[l - x] = N(ihTy m ^ 1 + x ^ = 2 nhn + i)2 1 - (3 - 30) 

For maximal transcedentality contributions this property also holds, in certain sense. It cannot 
be applied to integer N since the reflected image of, say, the leading harmonic sum S\(N) has 
poles at each positive integer point: 

Sl( -N-l) = Sl (N)-'-^- ("I) 

sin ttN 

However, the symmetry holds for half-integer values of N, from where the relation (|3.28p can 
be continued onto the entire complex TV-plane. 



3.4 Large x 

Quasi-elastic scattering, (1 — x) <C 1, is determined by the behaviour at large N, where 

S a (N) ~ N- a+1 , S- a (N) ~ \N~ a , (3.32a) 
y_ a (iV) = 0(N~ 2a+2 ) (a > 2). (3.32b) 

The second loop evolution kernel (|2.18bp is "regular" at N = +oo (that is, does not contain 
In N terms) and is dominated by the positive index term, 

S 3 (N)~N- 2 5V 2 (x) ~ (l-x). (3.33) 

It is down by two powers in (1 — x) counting as compared to the Born term a p h(l — x) _1 , while 
the negative index contribution is suppressed ever further. Power counting looks natural: the 
graph FigOJb) contains in the intermediate state two hard partons, each producing the phase 
space suppression factor d?k/k oc (1 — x) 2 , resulting in 

Y- 3 (N) oc N- 4 => 5V 2 (x) ~ (l-x) 3 oc (1 - x)~ l ■ [(1 - x) 2 } 2 . 

Unlike the second loop, the third loop contribution to the kernel V 3 exhibits a single power 
of logarithmic enhancement, Si(N) oc IniV. Now the positive index term is negligible, S5 oc 
./V -4 => a ph(l ~~ x ) 3 > an d ^ i s the last singular term in (12,18cP that turns out to be leading in 
this limit: 

iC2-S 1 (7V)S_ 2 (7V) ~ |^ 5p 3 (x) ~(l-x)ln(l-x). (3.34) 

It is worth noticing that the origin of the contribution (|3.34p is rather specific. Namely, it 
originates from harmonic sums with two negative indices. In QCD, maximal transcedentality 
sums of this nature in the third loop are [1] 

=> iC F (2C F -C A )[-2(C F +C A )S- 3 ,„ 2 -2(3C F -C A )S- 2 ,- 3 

+4C A< S_ 2 ,_2,i + 4(3C F -2C A )S- 2 ,i,-2 + 2{2,C A -2C F )S 1 ^ 2 ,-2\ ■ 

Equating the colour factors, C F = C A , one obtains a logarithmically enhanced contribution to 
the universal AT =4 SYM anomalous dimension, 

gjTuni =>" 5-3,-2 + 5-2,-3 ~ ^1,-2,-2 — 5-2,1,-2 — 5-2,-2,1 = — 2 5l(5_ 2 + ^i) • (3.35) 
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The term built of positive index sums, S1S4, is taken care of by the RRE, while the squared 
negative index sum, S^ 2 , propagates to the evolution kernel (|2. 18c|) : 



— ^SiS^ 2 — —^Si(S-2 — IC2) 2 — — fd ' Si ~ ^SiS-2 + h& SiS-2- (3.36) 

The first term on the r.h.s. participates in forming the a 3 correction to the physical coupling, 
the second is O (in N/N 4 ) . The last term is a cross-product of the subtracted sum, S- 2 oc N , 
and S_2(oo) = ~\(,2, and gives rise to a bizarre correction (|3.34p . 

The third loop contribution (|3,34[) looks as a radiative correction to the previous order 
kernel (|3.33p . Taken at face value, the presence of this logarithmically enhanced contribution, 
though subleading in the overall (1 — x) counting, negates our attempt to construct a perturba- 
tive scheme such that relatively large higher order corrections in all orders would be generated 
automatically. The fact that this proposal stumbled at the level of the third loop, calls for a 
deeper insight into the physical origin of contributions described by harmonic sums with two 
negative indices, in order to incorporate this eventuality into the "inheritance" framework by 
further elaborating the RR equation and its evolution kernel. 



3.5 Small x 

Small- x behaviour of the anomalous dimension, and of the evolution kernel, is determined by 
the rightmost singularity in the N plane. In reality, the BFKL limit in TV =4 SYM is described 
by the anomalous dimension 7 un i(A r — 2) having a pole at N = 1. For the sake of simplicity, 
we will ignore the shift of the argument and discuss the behaviour of V(N) at N = — 1 + lo for 
u — > 0. Participating harmonic sums have the following expansion at this point: 

Si(u) = -tv^ + b-u +^(-l) fc C^ fc ~\ (3.37a) 

k=3 

S a (iv) = -u>- a - Ca + 0{u) , a>2; (3.37b) 

S^ m (iv) = -S m (co) + O(u; ) , S- 2 (w) = u- 2 - fa + O(oj) ; (3.37c) 

Y_ m H = S- m (u;) -C 2 S- m+2 + ...; u; = N + l->0, (3.37d) 

where the factor (— 1)^ has been suppressed in S- m and Y- m . At small x, positive and 
negative index sums are equally important. 

According to (|3.37p . on the second line of (|2.18cp the leading and first subleading singu- 
larities cancel, 



F_ 4 -|(S_4 + ^ 2 ) + iC 2 5- 2 



- 1 c 2 ^ 2 ) - i( w - 4 + k 2 - for) + \t**- 2 = o{u-i) , 



(3.38) 



and the evolution kernel can be approximated as 



Vi = - Si ~ to' 1 - ( 2 uj + . . . , (3.39a) 

V 2 ~ ±S 3 - iy_ 3 ~ -to~ 3 + fou;- 1 + ..., (3.39b) 

V 3 -i55 + |y"_ 5 + i5 3 ~ 2oT 5 - 2C 2 w~ 3 + • • • • (3.39c) 
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To obtain the anomalous dimension one has to add to (|3.39p the generated pieces (jl.lOp . 
From (|3.39p one easily derives perturbative expansion coefficients for the generated reciprocity 



respecting (R) and violating terms (V): 

R 2 = 0, (3.40a) 

V 2 = ViPi = -uj- 3 + 0(1); (3.40b) 

Rs = ViVf + \V\Vx = 2uj- 5 - C 2 ^~ 3 + 0(1) , (3.40c) 

V 3 = ViV 2 + V1V2 = 4uj- 5 - 3( 2 uj-' 3 + 0(uj- 2 ) . (3.40d) 

The inherited contributions, 

R 2 + aV 2 = -auj~ 3 + 0(1) , (3.41a) 

Rs + vV 3 = 2(1 + 2a)uj' 5 - (1 + 3<t)( 2 lj- 3 + ... , (3.41b) 



also exhibit double logarithmic series and contain subleading singularities. Adding the genuine 
third loop evolution kernel (|3.39[) and the generated pieces (|3.41|) . for the space-like channel, 
a = — 1, we observe cancellation of double logarithmic terms, together with single logarithmic 
singularities: 

232 
7 (S) ^^ + .^ +0 .^i + C 2 ^ + ..., u, = iV + l«l, (3.42) 
oj uj z u) a 2lu 

the latter corresponding to two famous "zeroes" in the BFKL anomalous dimension [26]. 
For the time-like evolution, a = +1, combining (|3.39p and (13.411) produces 

7 (T) K ^_ 2 .^ + 8 .°k +C2 .^_ 6 ^|V.... (3.43) 

3.6 RRE and coherence in small-x physics 

Thus, in the small-x region the evolution kernel Vk contains double logarithmic enhancement 
factors, (a p h/w 2 ) fe ~ 1 , multiplying the one loop anomalous dimension, a p h/uj. In fact, the 
leading singular contributions given by the first terms in (|3.39p represent the start of the 
all-order series, 

r DL = ^(j^T^-J)^-^ + 2^-5^ + U% + ... i (3.44) 

Z \V r J UJ UJ 6 0J° UJ' UJ J 

that are characteristic for the evolution kernel. Its origin lies simply in the choice of the 
logarithmic ordering variable that separates successive parton splittings, known as "parton 
evolution time". As has been shown in [18], in order to unify the space- and time-like parton 
dynamics one has to choose parton fluctuation lifetime, k + /k 2 , as a common evolution variable 
for both channels. It is this choice that leads to the non-local evolution equation (jl.6p whose 
universal kernel V no longer coincides with the anomalous dimension. The latter acquires an 
additional term, o~V, which is opposite in sign in the two channels. This specific contribution, 
which is solely responsible for the breaking of the GL reciprocity, is of pure kinematical origin 
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and is easy to derive since, in a given order of perturbative expansion, it is inherited from 
lower orders. 

Inherited R and V contributions (jl.lOp also posses DL enhanced terms, see (|3.4ip . They 
combine with those of the kernel proper (|3.44p and cancel in the space-like case to produce 
single logarithmic BFKL series (|3.42p . The cancellation of DL enhanced terms in 7® holds in 
all orders, the reason being a general physical observation, due to V.N. Gribov, that inelastic 
diffraction vanishes in the forward direction. 

Indeed, consider a space-like scattering process in which parton p emits a soft gluon, p — > 
p' + k\ with k\+ <C p+, which gluon, in turn, radiates a still softer gluon k 2 (k 2 + *C The 
logarithmic condition, according to the fluctuation lifetime ordering, reads 

k 2 k 2 

h+ k 2+ 

This region includes, apart from k 2 4 < k| t , a ki t integration over the interval 

14 < kf t < |±i-k 2 4 , (3.45b) 

which is rather large in the soft kinematics and contributes as alns ■ a In 2 x to the evolution 
kernel V 2 . However, in the anomalous dimension this contribution cancels when one takes 
into account emission of k 2 off the partons p and p' . Physically, the process can be viewed 
upon as break-up (inelastic diffraction) of the incident particle p in the external gluon field 
of transverse size Aj_ ~ 1 / ^2* • However, in the kinematical region (|3.45bp the transverse size 
1/kit °f the parton fluctuation p — > p' + k\ is smaller than the resolution power of the probe, 
Aj_. In these circumstances the destructive interference between k 2 interacting with the initial 
(p) and final state (p 1 + k\) comes onto the stage and eliminates ()3.45b|) reducing lifetime 
ordering to the transverse momentum one, 

(V + R)-V : \4/k + =s> k 2 . 

In the language of the RRE, it is the role of generated DL terms to take care of this trans- 
formation due to coherent suppression of the part of the phase space (I3.45bp kinematically 
allowed to soft gluons by the dynamically unaware, unsuspecting lifetime ordering. 

For the time-like evolution, the very same generated double logs add to those in V, resulting 
in the characteristic DL series for the time-like anomalous dimension, 

7£? = 3 ( V^T^ - J) = 5* - 2^ + 8% - 40^ + 224^ + . . . , (3.46) 

whose first three terms are present in (|3.43p . This structure of the time-like anomalous dimen- 
sion is equivalent to ordering soft gluon cascades in angles: 

(V + R) + V : kf/fc + -> \4/k\. 

The angular ordering following, once again, from destructive soft gluon interference in the 
angle disordered kinematics [27], can be considered to be an image of the transverse momentum 
ordering in the space-like case, via the RRE. 
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(3.45a) 



In spite of the fact that the evolution kernel V does not correspond to a clever choice of 
the evolution variable in either T- or S- channel, its universality can be exploited to relate DIS 
and e + e~ anomalous dimensions. In particular, the RRE links together two puzzling results 
that were never thought to be of a common origin [28] . They are: the absence of the a 2 and a 3 
terms in the BFKL anomalous dimension (|3.42[) in the DIS problem and, on the other hand, 
phenomenon of exact angular ordering [29,30] that seems to hold, unexpectedly, down to the 
next-to-next-to- leading order in e + e~ anomalous dimension [31]. 



4 Conclusions 

Maximal transcedentality Euler-Zagier harmonic sums that describe anomalous dimensions 
of twist-two operators in TV = 4 super-symmetric Yang-Mills theory are genetically related, 
from QFT point of view, to multiple radiation of soft gluons. The word "soft" here stands 
for lack of better terminology. In fact, these gluons may be perfectly "hard", carrying a large 
(light-cone) momentum fraction x = 0(1). One would rather call them classical. The point is 
as follows. Mellin image of the Born level N '=4 SYM anomalous dimension, 

7i(iV) = -5i(7V)=V(A r + l)-V'(l) 710*0 = T~— — (z<l), (4-47) 

1 — x 

is nothing but the universal part of the gluon radiation probability which, according to the 
celebrated Low-Burnett-Kroll theorem [4] has classical nature. Such radiation does not depend 
on quantum properties of the emitter, but its (colour) charge, and does not affect quantum 
state of the radiating system. Classical (or LBK) gluons (14. 47ft constitute a significant part of 
diagonal QCD splitting functions, 



lg-*g{x) 



Cfcy 


X 




1 — X 


Caol 


X 


7T 


1 — X 



+ (1-x) 



2 



+ (1 - x) ■ (x + x" 1 ) 



(4.48a) 



(4.48b) 



The quantum parts of emission probabilities, depending on the details of the emitter, are 
down by two powers of the gluon momentum fraction (1 — x), in perfect accord with the LBK 
wisdom. 

"Classical" does not necessarily means "simple". However, amazing simplifications do 
arise, now and then, in the QFT context where classical radiation is concerned. Thus, exact 
Parke-Taylor QCD amplitudes [32] which describe radiation of arbitrary number of gluons 
with "classical" helicities, 2 — > n, and generalise production probabilities [33] of soft gluons 
with Xi <C 1 to arbitrary momentum fractions x ~ 1. Hence, a breathtaking possibility of 
an all- order prediction [34] of the magnitude of the "cusp anomalous dimension" [35], which 
is just another name for intensity of soft (classical) gluon radiation — the "physical" (or 
"bremsstrahlung" ) QCD coupling [2,3]. 

From this perspective, dynamics of the M = 4 SYM model has a good chance to be "classi- 
cal", to extent, since, thanks to super-symmetry, quantum effects due to fermions and bosons 
have cancelled both in the anomalous dimension (I4.47|) and in the /3-function. This idea finds a 
strong support in the conjectured gauge-string duality [36] — AdS/CFT correspondence [37]. 
A powerful property of integrability, both in the high energy Regge asymptotics of scattering 
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amplitudes [38, 39] and in the scale dependence of quasi-partonic operators [40] manifests it- 
self most profoundly in the Af = 4 SYM case [41]. Descending from SYM back to QCD by 
reducing the number of super-symmetries, one finds integrability property to survive in cer- 
tain sub-sectors of QCD parton dynamics, such as baryon wave function and maximal helicity 
multi-gluon operators [42]. 

All that makes this model a perfect playing ground for unveiling internal structure of 
multi-gluon radiation in QCD, whose effects spoil perturbative expansions and bear major 
responsibility for complexity of higher order calculations, and of QCD anomalous dimensions 
in particular. In the present work we chose J\f = 4 SYM model in order to learn how to 
single out higher loop effects, and structures, that are inherited from the lower orders of 
the perturbative expansion. To this end we used bookkeeping based on the concept of the 
Gribov-Lipatov reciprocity respecting evolution equation, RRE, [2, 18] which treats space- 
and time-like parton multiplication in one go. 

Power of the LBK wisdom about the classical nature of radiation described by (|4.47|) re- 
vealed itself most remarkably in this model. The RRE applied to the space-like anomalous 
dimension 7 un i, known in three loops [5], has generated the major part of higher order con- 
tributions in terms of a compact "evolution kernel" V linked to 7 un ; by a non-linear relation 
(jl.8p . Simplicity of the evolution kernel is suggestive and invites deeper studies. 



Note added. When this project was under completion, we learned about the work by 
B. Basso and G. Korchemsky in which the reciprocity respecting equation (jl.8p has naturally 
emerged as a consequence of the conformal symmerty. 
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A Harmonic sum representation 

The answer contains harmonic sums, 

N N 

S a = Yl k ~ a > S - = S P (N) = S P (N) - S p (oo), 

k=l k=l 

having the following Mellin representations: 

x 



Si (AT) = M 



1>{N + 1) - V(l) , 



and (a > 2) 



Sa(N)-Ca= S a (N) 

S- a (N) + ( a (l-2- a ) = S- a (N) 



(-1) 



a-l 



r(o) 
r(o) 



M 
M 



x ln a 1 x 
x-1 

xln a_1 x 
x + 1 



(A.la) 

(A.lb) 
(A.lc) 
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As for generalized multi-index Euler-Zagier harmonic sums, a representation alternative to 
that in f)2 . 13|) is better suited for analysis of the large- N asymptotics. It employs the basis 
of sums with the negative index moved from the head to the tail of the index vector, 5 r ^ j _ p 
instead of S_ p ^. The latter may contain a logarithmic enhancement factor, 5_ p ,,^ oc S\{N), 
with £ the number of the rightmost +1 indices, e.g., 5_ 2,1,1,1 SfS-2 = 0(ln 3 iV/iV 2 ). 

Crucial simplification of the answer is achieved in terms of complementary harmonic sums 
Sm-p which we define by a recursive relation 



., an (N) = S ai _ an (N) 



n-1 

^"1 Sai,...,a k 
k=l 



(AO • 5, 



,a„(oo), n > 2. 



(A.2) 



In particular, 

Sg,- P (N) = S a ^ b (N) - S a (N)S- p (oo), 
Sg, b ,- P ( N ) = S a ,b,-p(N) - 5 a (AQ5 6 ,- p (oo) - 5 a , 6 (AT)5- p (oo), 
S a ,b,c,-p( N ) = S a,b,c-p( N ) ~ Sa(N)S bjC - p (oo) - S atb (N)S c - p (oo) - S ajbtC (N)S - p (oo). 



These functions are finite at iV = 00. Subtracting the value at infinity, one arrives at (|2.16p . 
The case of all positive indices aj but the very last one, a n = —p, is special in that such sums 
fall fast with N and contain no lnA^ factors accompanying any subleading power N~ k . 
Mellin transforms of participating multi- index sums are as follows: 



rs, 



TS a . 



a,—p 



b,—p 



rs, 



a,b,c,—p 



where 



M 
M 
M 



X 


r dz 1 


1 + x 


Jx z + 1 


X 


r 1 dz 


1 + x 


J x l + z 


X 


r 1 dy m 


1 + X , 


L 1+1/ 



In 



a-l 



ln^ 1 - 



x 

In?- 1 - 

z 

a-l y 



Z 

dt 



In 6 " 1 - In"" 1 - 



1 + t 



a-l 



t 



(A.3a) 



dt 



X Jy 1+t 



In 6 - 1 - 



, ^l^" 1 - ln^ 1 - 

y Jt l +z t z 



-1 



n ^ 

' se(a,...,p) 

Constructing linear combinations (|2.14dp using (|A.3ap . one arrives at ()2.14cp . 



(A.3b) 



B Shifting negative indices to the right 

Here we list relations that help to transform the original representation (|2.13p in terms of 
harmonic sums S- p ^ of [5] into the negative-on-the-right form, S^-p'- 
Two indices 

5—2,1 = — Si,— 2 + S\S-2 + 5— 3) S-2,2 = ~ $2,-2 + S2S-2 + S_ 4, 
S-2,3 = —S3 t -2 + SsS-2 + S-5] S-3 t i = —Si } -3 + S1S-3 + 
5-3,2 = —52,-3 + 525-3 + 5-5; 5_4,i = — 5l,_4 + 5i5_4 + 5_5. 
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Three indices 



5-2,1,1 
5-3,1,1 
-5-2,2,1 



+ 



Si, 1,-2 + Si(S-3 — Si^2) + 5'_25'l ) l — 5*2,-2 — S\ -3 + S_4, 
Si, 1,-3 + 5'i(S'_4 — Si^s) + S-^Si^i — S*2,-3 — 5^-4 + 5_5; 
5-2,1,2 = (52,1,-2 + -51,2,-2) + 5'i(S'25_2 — £2, -2 + ^-i) 

+ S21S-3 — Si : -2) + S3S-2 — 2S , 3,_2 — -52,-3 — 5l,-4 + 25'_5. 



Four indices 

3-5-2,1,1,1 



— — 35i,i,i,-2 + 5i,i(25_3 — 3«5i,_2) + 5i,i<5_3 + 5s5_2 
+ 51(3)51,1,-2 + (5i,i + S < 2)5_2 — 3(52,-2 ~~ 35i,_ 3 + 3S'_4) 

+ 3(<5l 2,-2 + 52,1,-2 + 51,1,-3) — 3(j5i -4 + 52,-3 + 53,-2 — 5_5) 



[hep-ph/0403192] 



[hep-th/0404092] 
|hep-th/05020T5] 
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